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1. Introduction 

The virial theorem has a wide range of applications in physics. For example, the virial 
theorem is a powerful tool for studying static, non-evolving (relaxed) systems such as 
stars, gas clouds, star clusters, galaxies and galaxy clusters. It can also be used to study 
large structures in the universe. Not only is it applicable to nonrelativistic dynamical 
systems, but it can also be formulated to deal with relativistic systems which require 
special relativity or even general relativity for their description. Of course, one cannot 
expect to obtain a complete description of the system only by using the virial theorem. 
However, this theorem can provide considerable insight into the dynamical behavior 
of the system. For example, Zwicky could discover the large mass discrepancy in the 
galaxy clusters only by applying the virial theorem to the Coma cluster [TJ. 

For a historical review of the virial theorem and also its applications in stellar 
astrophysics, see [2]. 

In the context of the alternative metric theories of gravity, it is known that only the 
values of a set of numerical coefficients in the components of the post-Newtonian metric 
vary from theory to theory. Thus, one can encompass a large class of alternative theories 
by introducing arbitrary parameters in place of the numerical coefficients. In the other 
words, this framework provides a model independent tool for testing the metric theories 
of gravity in the solar system in the small velocity and weak-field regime. Although 
this idea dates back to Eddington [3], but the parametrized post-Newtonian (PPN) 
formalism was developed by Nordtvedt and by Will (see jl] for more details). 

In this paper, we derive the PPN virial theorem. In fact, our main result is a 
post-Newtonian virial theorem which can be used in a large class of metric theories of 
gravity. It should be noted that the post-Newtonian virial theorem in general relativity 
has already been investigated by Chandrasekhar [5j [6] . 

The outline of this paper is as follows: In section [2} we derive the PPN 
hydrodynamics equations for metric theories of gravity. In section^ we derive the tensor 
and scalar form of the PPN virial theorem from the PPN hydrodynamics equations. 

2. Hydrodynamics equations in the post-Newtonian approximation 

Einstein equivalence principle is the foundation of metric theories of gravity. Thus 
in these theories, the world lines of test bodies are geodesies of the metric g^ v . 
Mathematically, this postulate can be written as the following identity 



where T^ v is the energy-momentum tensor. We assume that for a hydrodynamic system 
the energy-momentum tensor is given by 



where p is the rest-mass-energy density of atoms in the fluid element, p is the pressure, II 
is the specific density of internal kinetic and thermal energy in the fluid element and 



V^ v = 
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is the four velocity of the fluid element in the comoving frame (we use the geometrized 
units in which the vacuum speed of light and the gravitational constant are unity i.e. 



c = 1, G = 1). In the post-Newtonian limit, 
gravitational potential) are small quantities: 



dt > 



v 2 ~ 



u~ n 



£ ~ 0(2) 
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where 0(n) stands for a term of order i where c is the vacuum speed of light. On the 
other hand the PPN metric is given by [I], 

g m = -1 + 2U - 2/3U 2 - 2^ w + (2 7 + 2 + a 3 + (i - 2f)$i (4) 
+2(3 7 - 2/3 + 1 + C 2 + 0^2 + 2(1 + Cs)$3 
+2(3 7 + 3C4-20$4-(Ci-20^ 



#0; 



3 + «i — a 2 + Ci 
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^ = (1 + 2 7 C/)% 

where f7 is the Newtonian gravitational potential determined in terms of p as 
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and ^l) $2, ^3, $4, Vj, Wj, and *4 are other metric potentials (see Appendix A 
for the definition of these potentials). a±, a>2, «3, Ci; C2, C3> C45 C> 7 an d are the PPN 
parameters. As it is known, each parameter has a special meaning and one can compare 
different metric theories in the post-Newtonian approximation only by comparing their 
PPN parameters. 

Energy-momentum tensor T^ v and the Christoffel symbols T^ u can be written to 
the post-Newtonian order by using the PPN metric. See table 4.1 of pi] for in the 
post-Newtonian limit and table 6.1 for Christoffel symbols. Now, by substituting the 
metric components (SI)- (EJ) into equation (TI|) and by using the tables 4.1 and 6.1 of [I], 
the time component of (CQ) takes the following form 

da d(av 



dt 



dx i 



+ (3 7 - 2 ),f 



dp , .dU i 
^ + (37-3)^ 







where 
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p(l + U + v 2 + 2U + 



p' 



(8) 



(9) 



and v = is the velocity of a fluid element. Equation flHj) is the continuity equation 
in the post-Newtonian approximation. It is straightforward to show that the space 
component of (CQ) takes the form 



djcrvj) 
dt 
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where 



(<7 + fnv 2 )Q-- Ft [c 1 w i + c 2 v i \ = o 



c 2 = -(47 + 4 + ai - a 2 + Ci - 20 



Ci = -(l + « 2 -Ci + 20 



(11) 
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and 



$ = -(2 7 + 2 + a 3 + Ci - 20^i + (1 + Cs)$3 

+ (3 7 - 2/3 + 1 + C 2 + 0*2 + (37 + 3( 4 - 20*4 
In order to simplify equation (jTUJ) let us introduce the potential x as (01) 

-y p(x')ix-x'irfv 



X -- 

one can easily verify that 

d 2 X 



V s - W j 



(13) 
(14) 

(15) 



dtdxi " J ' 3 

on the other hand, for a perfect nonviscous fluid the Euler's equation in the Newtonian 
limit is given by 



dVi dU dp 



(16) 



where 0(4) is the post-Newtonian corrections to the Newtonian Euler's equation. Note 
that we do not need the exact form of the post- Newtonian corrections 0(4). Now, by 
using equations (|T4l) - (|T6|) . equation (JTUJ) takes the following form 

"V + + + 3(7 - 1)u)p] - P M 



+p 

4 







I Qtn Cl 1 
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dx 

4 7 + 4 + a x .dVj 
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(17) 



+(5 7 -l)p 



eft 



v,U 
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<9af 2 cte* 



where 



X 1 + 7 2 1 n , 3 7 -2/? + l rr 3 7 p 

= it H U H 

^ 2 2 2 2 p 

V 2 <f> = -87TP0 



(19) 



Equation (|17p is the Eulerian equation in the post-Newtonian approximation. Equations 
(117)1 and (jSJ) are enough for describing a hydrodynamic system in the post-Newtonian 
approximation. 
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3. PPN virial theorem 



In Newton's theory of gravity, the total energy of a gravitating system is conserved. 
Furthermore, if the gravitating system is not affected by external forces, its total 
linear and angular momentum are also conserved. However, in the post-Newtonian 
approximation this is not the case for all metric theories of gravity pE]. For example, for 
the conservation of the total linear momentum of system, the parameters Ci>C2)C3>C4 
and «3 should be zero [4]. It is obvious that, for theories which violate the conservation 
of the total linear momentum, we can not find a PPN virial theorem. Because in these 
theories, the given gravitating system would not be stationary or periodic in time in the 
post-Newtonian approximation and consequently it can not be considered as a virialized 
system. Thus we shall restrict ourselves to metric theories in which five PPN parameters 
Ci, C2, Cs, C4 and a 3 are zero. 

Now, let us rewrite equation ffTT|) as follows 



djavj) 
dt 



+ 



8 . 



+ Jrf(l + 3( 7 -l )E 0p]-*g 



(20) 



1 + ot 2 
+^ — P 



Wi) - W ? ; 



+ 



47 + 4 + ct! jdVj 



dx l 



+(5 7 -l)p 



d 
dt 



47 + 4 + 01 " 
ViU — —Vi 







2(5 7 -l) 

A - $ 2 + 2$ 4 

where we have used the following expressions 

B = x,oo -A + $1 




jiVi-Wi] 



(21) 
(22) 



In order to find the PPN virial theorem, we multiply equation f l20|) by Xj and integrate 
over the volume V of the fluid. Most integrals can be calculated by using the lemmas 
introduced in [6]. So we do not present the detail of the calculations here. However, 
there are four new integrals which we combine them in the following expression 

d 



Mi. 



- I x 3P-q^[ b + ®w - $2 + 2$ 4 ]d 6 x 



(23) 



in the Appendix A we have shown that My is a symmetric tensor. Hence, the result 
can be written as 



dt 1 3 



2%j + m {j + s^p - 



47 + 3 + a\ - a 2 



(24) 
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and J Uid 3 x is the conserved total linear momentum of the system [4]. In fact, by 
integrating equation (120]) over the volume of the fluid, one can verify that 4 J Hid 3 x = 0. 
The tensors on the right-hand side of equation ( 1241) are defined as follows 
1 



Tij = - J OViVjSx 

5 lj P = S il J [(l + (3^-l)U)p]d 3 x 
dU 



3 dx l 



Mi 



pXjV 



X jP 



k ,3 



dx i 



d x 



■ dU 15$ 

dx % 2 dx l 



d 3 x 



Qij 



pUviVjd 3 x 
pVjVid 3 x 

p(x)p(x')v(x').(x- x 





p(x)p(x')v(x).(x — x')v(x').(x — X 

3Cx-x'Ux-x" 



/N ^ (x)(a: j - a^j + Vj(x)(xi - ^ 33, 

I v — y' I 3 



A [ x i x i)( x j x j) ,3 ,3 / 

—dxdx 



r'\5 



(26) 
(27) 
(28) 
(29) 

(30) 

(31) 
(32) 
(33) 
(34) 



M l3 = -j J J p(x)p(x>(x")[(^- 



fx-x"Ux'-x" 



|x — x'| 3 |x' — X' 



(rp . rp ^ ( rp . rp \ I ( rp . rp ^ ( rp . rp^ ^ 

7 ! - ; - ' ii\nf' , .r'(i'\r" 



x — x 



/|3| Y ' 



(35) 



These tensors are symmetric in i and j (see [B] and Appendix A for more details). Since 
the right-hand side of equation ff24|) is symmetric then the antisymmetric part of the 
left-hand side of (J23J should vanish, i.e. 

d 



dt 



J (xjUi — XiUj)d 3 x = 



(36) 



this equation expresses the conservation of the total angular momentum of the system. 
If the system is stationary or periodic in time then the time average over a long time 
interval of the symmetric part of the left-hand side of equation (1241) vanishes, i.e. 



j t J(x j n i + x i U j )d 3 x^ ~ 



(37) 



where (...) denotes the time average. Thus the PPN generalization of the classical tensor 
virial theorem takes the following form 

2 {%) + (0%) + * y (P) - 47 + 3+ o ai ~" 2 (u* + %} + 2 + (5 7 - 1) (Wij) (38) 
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By contracting i and j, the scalar PPN virial theorem is 
2(T) + (fflJ) + 3(P)-ll±i±^i(u) 

+2 (4) + (57 - 1) (W) - (2) + £ (M) = (39) 

where 

(T> = \(J wWx^ (40) 

W = -\(J pUd 3 x^ (41) 



(W) = <^J pv l Ud 6 xj (42) 
(u) = -2 (it) = ( / ptv^ 3 x\ (43) 



<f>^ = - (^p0Ud 3 x^ (44) 
«) = 2 <Z> = 2 (/ / p(x),(x') VW - (X - x X l V ff- (X - X,) ^) (45) 

(M> = -2 (y 1 1 p(x)p(x>(x") ( * " ^ " d?xM^ (46) 

Equations fl38|) and fl39|) are our final results. It is clear from equation (1391) that the 
scalar PPN virial theorem has six extra terms relative to the Newtonian virial theorem. 
If we set 7 = /3 = 1 and ol\ = «2 = £ = 0, then we obtain the post-Newtonian virial 
theorem in general relativity. In this case, our result is identical to the result obtained 
by Chandrasekhar [6]. 

It is obvious that the PPN virial theorem can be applied to investigate the post- 
Newtonian effects of metric theories of gravity on the hydrodynamic behavior of the 
systems which are in the post-Newtonian limit. By applying the PPN virial theorem to 
such cases, one can compare metric theories of gravity with general relativity. In these 
systems, general relativistic effects on the dynamics cannot be ignored. For example, 
post-Newtonian corrections in supermassive stars and massive white dwarfs would lead 
to unexpected gravitational instabilities. This fact has been shown by means of rather 
detailed calculations [7J [8] . However, by using the post- Newtonian virial theorem this 
result could be anticipated without the need of detailed calculations [21 [9]. Of course, 
as we mentioned before, by using the virial theorem one cannot expect to obtain as 
complete a description of a hydrodynamic system as would be possible from the solution 
of the field equations. However, the virial theorem can provide extensive insight into 
the behavior of the system. 
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4. Conclusion 



In this paper, we have derived the PPN virial theorem. The virial theorem in the post- 
Newtonian approximation has already been derived only in general relativity. We have 
generalized it to every metric theory of gravity. Albeit, we mentioned that the PPN 
virial theorem can only be written for theories in which parameters £1, (2, C3> C4 an d 0^3 
are zero. 



Appendix A. Tensor Mj. 



In this appendix we simplify My (integral ( 123]) ) as far as possible. Also we show that 
My is a symmetric tensor. First, let us write the metric potentials of the PPN metric 
(see [1] for more details) 



$1 = 

$ 2 = 
$ 3 = 
$4 = 

A = 
B = 



x — X 



p'p" 



x' - x" 



x — X 
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x — x' 


3 ' 




x- x" 




x' - x" 





d Xj d oc 



|X — X' 
JTTI 



SSL** 



|x — X' 

p'W 
|x — x' 



V 



x — X' 



-d 3 x' 



-d 3 x' 



p'[v'.(x-x')] 2 



x — X 



t\3 



-d 3 x' 



|x — x'| dt 

P' V 'j j3 f 

J —d 4 x 



x — X' 



pV.(x-x / )(x 3 --^) 

-d x 



x — X 



'13 



(A.l) 
(A.2) 
(A.3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 



where prim means /' = /(x'). 

In order to simplify equation (I23I) . we start with its first term. Using the definition 
of B and the Euler's equation (JIB]) we get 

Xjp dp 



dB, 3 

X n O ~ . Ct X 

3 dx l 



Xj(xi - a?i),_ ^ M dp 



|x — x'l dx'' 



|x — x i 



,, 3 ( x k x k) 



dx' k 



d Xj d x 



+ pp'p' 



1 ji x j ( x i x i)( x x )-( x x 1 d z xd 3 x'd z x" 
|x — x'l |x' — x" '" 



/ // 

ppp - 



A 1 i] d 3 xd 3 x'd 3 x" 



x — x' 1 1 x' — x" 1 3 



(A.10) 
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using integration by parts, one can easily show that the first integral on the right-hand 
side can be written as 

p(x') d ( Xjp(x') 



2 / p 



Si 



x — X' 



dx i 



X — X 



(All) 



on the other hand, by using the definition of $4, the above equation is equal to 

~ 2 f X ^ d3x (A - 12) 



Thus we can write 

f 96 J3 o f 



d® 4 

dx % 



+ PPP 



- / PP'P"- 



'~" Xj ( Xi ^ )(X X/) ' (X ' *"±d 3 xd 3 x'd 3 x" 



|x — x'l |x' — x"| 3 



(A.13) 



ry . I ry' ry | 

X — x' 1 1 x' — x" 1 3 



d 3 xd 3 x'd 3 x" 



Now consider the second term of equation 
straightforward to show that 

,3 



Using the definition of $vk it is 



y ^iP^T" 6 ^ = / Pp'p"x j H i (^x.' 1 x.")d 3 xd 3 x'd 3 x" + J Xjp^fd 3 x (A.14) 



Z?)(x X').(X' X"]_ d 3 xd 3 x , d 3 x/I 



x — x'| 3 |x — X' 



//|3 



where f/j(x, x',x") is given by 



3(x< -xO(x-x') /(x-x") (x'-x"; 



|x — X 

+ 



r\5 



{A - O 



r'|3| 



|X — X" 

|x — x'| 3 |x' — x"| 



x x x x < A ' 15) 

it is obvious that Hj(x, x',x") = — f/j(x', x, x"). Now we can simplify the first integral 
on the right-hand side of (IA.14I) . Since the integrand is symmetric in x and x' (and also 
in x" and x') then the integral can be rewritten as 



I PPP 

PP'P 



, ff 3<-(x i -x / i )(x-x' / ).(x / -x" 



|x — x'Plx' — x"| 



-d 3 xd 3 x'd 3 x" 



(A.16) 



'»" [ 3 ] ' i \ d 3 xd 3 x'd 3 x" 



|X — X'| °|X — X' i 

Similarly the third term on the right-hand side of flA.14j) can be written as 



PP P 



xi)(x x").(x' j r , if r „ 

x — x' r x' — x" 3 



Also, the second integral on the right-hand side of (1A.13j) can be replaced by 

rpM ( ry ry \ 

J Ji 3 ^ ' *' 



PPP 



|x — x'| 3 |x' 



:d 3 xd 3 x'd 3 x" 



(A.17) 



(A.18) 
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Now, by using equations (1A. 13j) -( 1A.18[) and (1251) we get 

r f f , „3(xj — x'Mxj — x'-)(x — x').(x — x") „ „ , „ „ 

Ma = - / / pp'p"— ^rr tt -d 3 xd 3 x'd 3 x" A.19 

III ■ x — x" x — x' |b 




/ ^iX^j x j)( x x ).(x x ) „ „ „ 
pp p j 75— —5 a xd x ax 




|x — x'| J |x' — x , 

/ //(^ — ~ ^j) ( X i ~~ x j)( x i ~ x i)-, ,3 ,3 , ,3 „ 

po p : — —. p \d xd x d x 

'x — x'| d |x' — x"| 

it is obvious from this equation that is a symmetric tensor. 
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